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Gluodynamics and two-flavor QCD at non-zero temperature are studied with the so-called over- 
improved cooling technique under which caloron solutions may remain stable. We consider topo- 
logical configurations either at the first occuring stable plateau of topological charge or at the first 
(anti)selfdual plateau and find the corresponding topological susceptibility at various temperatures 
on both sides of the thermal transition or crossover. In pure gluodynamics the topological suscep- 
tibility drops sharply at the deconfinement temperature while in full QCD it decreases smoothly at 
temperatures above the pseudocritical one. The results are close to those calculated by other meth- 
ods. We interpret our findings in terms of the (in)stability of calorons with non-trivial holonomy 
and their dyon constituents against overimproved cooling. 
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I. INTRODUCTION 

Topological charge has many faces. It is well known 
that the appearance of its space-time distribution 
strongly depends on the method by which it is stud- 
ied. If resolved by means of the modes of a lattice 
Dirac operator with maximal chiral symmetry (the 
overlap Dirac operator), the resolution can be dialled 
by the cutoff applied to the eigenvalues of the modes 
(in a symmetric band enclosing the zero eigenvalue) 
included in that analysis. In this way very differ- 
ent types of topological structures are revealed, rang- 
ing from globally extended, laminar sheets of alter- 
nating sign to instanton-like lumps appearing in the 
infrared (IR) pHa]. To some extent the scale of ul- 
traviolet filtering can be mimicked by the number, 
say, of overimproved stout link smearing steps |^|. 
Before the fermionic methods became popular, only 
cooling Q and smearing pj were available for investi- 
gating the topological vacuum structure. This search 
was biased in favor of detecting instantons. Results 
concerning the "instanton structure" have been pre- 
sented and discussed in the "Confinement and Topol- 



ogy" sessions at the annual Lattice conferences un- 
til 2000 H-EH. Beginning from 1998, for non-zero 
temperature calorons with non-trivial holonomy and 
their dyon constituents [Il| - [l5| have attracted more 
and more interest, although within a relatively small 
community 16-29]. 

There was always the hope to elucidate the dif- 
ferent phases, both of pure Yang- Mills theory and 
of full QCD, in terms of the topological structure. 
The local behavior of the laminar sheets (and the 
resulting two-point function seen with few cooling 
steps [J, |30j), however turned out to be not criti- 
cally dependent on the phase. Only in the IR, at 
best, the characteristic differences may become vis- 
ible [33, [32[. One prominent example is the space- 
time anisotropy of the susceptibility in slab-like sub- 
volumes [33J] originally predicted in the instanton- 
antiinstanton molecule model [341 ] . 

In this paper we will use a specific (overimproved) 
kind of cooling (22|, [35|, however in the Cabibbo- 
Marinari mode for SU(3), and we are going to use 
it far beyond the point where the IR structure (as a 
mixture of instantons and antiinstantons) usually has 



been studied. Actually, this kind of action (see Eq. ([2]) 
below) was invented J_35] exactly for stabilizing instan- 
tons or sphalerons [35, 36] under cooling. A similar 
purpose has been pursued in Ref. |22j in order to study 
the elusive instanton constituents at T = 0. 

We will concentrate on the non-zero temperature 
case. We employ overimproved cooling in order to 
see in as far - after eliminating all short-range fluc- 
tuations - the emerging topological objects ("multi- 
caloron" configurations) will clearly distinguish be- 
tween the different phases or thermodynamic states 
in the neighbourhood of the deconfinement transition 
for gluodynamics and crossover for two-flavor QCD, 
respectively. The nature of these topological config- 
urations in the case of SU(3) gauge theory has been 
carefully considered already in Ref. [24| . and earlier 
for SU(2) gauge theory in Ref. [la ]. 

In this stadium of cooling the fields are either self- 
dual or antiselfdual or trivial throughout the lattice, 
i.e. cooling has been employed until it reaches the 
scale set by the whole lattice. Hence no more details 
than the total topological charge Q can be read off 
and be associated to the original, thermalized config- 
uration. This identification is convincing, however, 
since the topological charge mostly stabilizes after a 
few cooling steps. No further localization [32, |37| of 
topological charge characterizing the original Monte 
Carlo configurations is possible in the final stadium of 
cooling. 

Our study instead focuses on the following ques- 
tions: Under which circumstances nontrivial topolog- 
ical structures are stabilized with respect to cooling ? 
How does this depend on the confining / deconfining 
nature of the ensemble the gauge field configurations 
are taken from ? What is the influence of the average 
Polyakov loop on the result of cooling ? 

Our previous caloron studies [16l |24J , starting from 
thermalized (Monte Carlo) lattices at strong coupling 
or deep in the confinement phase, did not consider 
the role of temperature in detail. Next, at an inter- 
mediate scale of resolution, the study of topological 
objects in SU( 2) lattice fields at non-zero tempera- 
ture [25, 26, 28] has shown the profound difference be- 
tween the confinement and deconfinement phase. To 
be concrete, in this analysis either 50 APE smearing 
steps with a smearing parameter 0.45 [25] or an over- 
lap fermion analysis based on 20 lowest modes |2(| [28| 
have been applied. 

The following picture of the topological content of 
SU(2) lattice gauge theory has emerged. At low tem- 
peratures one finds topological objects represented 
by non-dissociated calorons with maximally nontrivial 
holonomy [ll| - |l3j |. With increasing temperature [38J 
their composite nature in the form of monopoles be- 
comes recognizable. They start to dissociate into 
dyons of topological charge ±1/2 that appear (in the 
limiting case) as static U(l) monopoles in the maxi- 
mally Abelian gauge. Approaching the transition tem- 



perature T c (to the deconfining phase) from below, 
approximately one half of the calorons were observed 
dissociated, retaining the symmetry between the con- 
stituent dyons. 

Above the transition temperature, a non-zero ex- 
pectation value of the averaged Polyakov loop devel- 
ops, which induces an asymmetry between the con- 
stituent dyons. This can be comprehensively ex- 
plained in terms of the peak values of the local 
Polyakov loop [28( : "light" (anti)dyons, with the lo- 
cal Polyakov loop of same sign as the averaged (pre- 
vailing background) Polyakov loop, become the most 
abundant topological objects, while "heavy" dyons or 
antidyons with the local Pol yak ov loop opposite to the 
average are suppressed J3JJ [4y] . The former "heavy" 
dyons may still carry highly localized fermionic zero 
modes bound to the "defects" detectable in the local 
Polyakov loop (4J, |42| • Nondissociated calorons are 
even more suppressed. 

First steps towards the statistical mechanics of self- 
dual dyons for the case of SU{2) gauge theory, sup- 
posed to be valid in the region around the critical tem- 
perature and derived from the semiclassical partition 
function, have been made in the papers [43l |44| . 

The present paper does not attempt to consider 
topology at the above mentioned intermediate scale 
of resolution. Concerning SU(3) gauge theory, this is 
left to an overlap fermion analysis of gauge field con- 
figurations which is in progress. The main outcome 
of the present study will be that the phase transi- 
tion between confinement and deconfinement in SU(3) 
gluodynamics can also be characterized by the sharp 
change in the appearance of topological objects re- 
maining once cooling has hit the lattice scale. We 
confront this with corresponding (softer) results for 
the crossover known to replace the phase transition in 
the case of full QCD with Nf — 2 dynamical fermions. 

We claim that our observations concerning the be- 
havior of action and Polyakov loop can be explained 
by the dyonic structure of calorons. We investigate 
the volume and discretization effects for the topolog- 
ical susceptibility calculated by the use of topologi- 
cal charges measured either at plateaus of topolog- 
ical charge or identified by the coincidence between 
\Q\ and S/S- Uls t and compare it with the topological 
susceptibility of uncoolcd Monte Carlo configurations 
calculated by other methods [45M471 ] . 

The following Section [TT] contains the main defini- 
tions as well as some details of the simulations we 
have used here or where the full QCD configurations 
are taken from. Section IIIII shows examples of cool- 
ing histories for SU(3) gluodynamics in the confined 
and deconfined phases. We comment on the influence 
of the overimprovement parameter. The temperature 
dependence of the calculated topological susceptibil- 
ity, both for pure gauge theory and for QCD, is dis- 
cussed in Section lTVl In Section|V]the evolution of the 
averaged Polyakov loop during the cooling process for 



SU(3) gluodynamics and for full QCD, both in the 
confined and deconfincd phases, is presented and the 
interpretation with the help of the dyonic picture of 
the gauge field ensemble is discussed. Section |Vl] is 
reserved for conclusions and discussion. In the Ap- 
pendix we recall some facts on SU(3) calorons and 
dyons. 



II. THERMAL ENSEMBLES 

For the study of the phase transition in the SU(3) 
pure gauge theory we employ the standard Wilson ac- 
tion Sw with the lattice coupling j3 — 6/ '<?q where go 
is the bare coupling constant. To determine the corre- 
sponding lattice spacing a as a function of /3, for this 
action we have used the Necco-Sommer parametriza- 
tion [48[. In what follows we will refer to this case as 
gluodynamics. 

To study the topological aspects of the phase transi- 
tion or crossover in a theory with dynamical quarks we 
have studied gauge held configurations generated with 
the gauge action Sw and Nf — 2 dynamical flavors 
of nonpcrturbatively 0(a) improved Wilson fcrmions 
(clover fermions). The configurations had been pro- 
duced by the DIK collaboration [49| using the Berlin 
QCD code (BQCD) [50|]. The improvement coefficient 
csw was determined nonperturbatively [51|. The lat- 
tice spacing and pion mass has been determined by 
interpolation of T = results obtained by QCDSF 
[52j ] . In the remainder of the paper this case will be 
referred to as full QCD. 

Our calculations were performed on asymmetric lat- 
tices with the four-dimensional volume V = a Lt ■ L s , 
where L t is the number of sites in the time (4th) direc- 
tion. The temperature T is given by T — l/aL t . In 
the case of gluodynamics we have employed L t = 4, 
L s = 16 as well as L s — 24 lattices, for which we 
have generated and analyzed 1000 and 500 configura- 
tions, respectively, at a set of /3 values. The coupling 
fi c = 5.692 [53[ characterizes the transition at L t = 4. 
It cor resp onds to the critical temperature T c ~ 300 
MeV |46j |. In order to study finite lattice spacing ef- 
fects we have also simulated L t = 6, L s = 24 lat- 
tices with statistics of 500 gauge field configurations 
for each /3 value. The phase transition in this case 
takes place at j3 c = 5.894 [53]. In order to keep au- 
tocorrelations small, all measurements were made on 
configurations separated by 500 sweeps. 

In the case of full QCD [49( we have analysed config- 
urations produced on lattices with L t = 8 and spatial 
sizes L s — 16 (500 configurations at each T/T c ) and 
L s = 24 (200 configurations at each T/T c ). The tem- 
perature T was effectively varied at fixed /3-value by 
changing the Wilson fermion hopping parameter n, 
i.e. the quark mass was not kept constant. The chiral 
crossover temperature T c sa 230 MeV was determined 



in Refs. |49|, [5J] at a pion mass value of 0(1 GeV). 



III. COOLING HISTORIES 

Technically, we cool down each gauge field configu- 
ration by means of the usual Cabibbo-Marinari cool- 
ing procedure, now with respect to an overimproved 
action [35j that has been used before [24J for the study 
of SU(3) calorons and multicalorons. We monitor the 
cooling process and search for plateaus of the action 
and topological charge appearing in the cooling his- 
tory. We consider 

• either the first plateau of the topological charge, 
when the topological charge Q stays near some 
integer value n (\Q — n\ < 0.1) for at least 100 
cooling sweeps, 

• or the first (anti)selfdual plateau, when the dif- 
ference between \Q\ and the action S (measured 
in units of instanton action Si ns t) becomes small 
(S/S inst - \Q\ < 0.1). 

For both gluodynamics and full QCD in the con- 
fined phase we observe the two kinds of plateaus to 
coincide, whereas deep in the deconfined phase no 
plateaus of any kind are seen. In the transition re- 
gion (just above the transition to the deconfined phase 
or for full QCD above the crossover) the first topo- 
logical charge plateaus can be different from the first 
(anti)selfdual plateaus. Moreover, both plateaus are 
unstable. 

We use overimproved cooling [35j because it par- 
tially stabilizes calorons with respect to shrinking and 
falling "through the meshes of the lattice" . The fol- 
lowing parametrization of the lattice action ( *-+w = 

U(t(x) S SU(N) link variable) allows for over- and 
underimprovement 
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Expanding in powers of the lattice spacing a one 
finds 35], 
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(note that no implicit summation convention is im- 
plied in this formula). S(e — 1) corresponds to the 
Wilson action, and the sign of the leading lattice ar- 
tifacts is simply reversed by changing the sign of e. 



Based on a discretized continuum one-instanton solu- 
tion of size p, one finds 
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suggesting that under cooling p will decrease for e > 
and increase for e < 0. 

In order to illustrate the influence of the parameter 
e, we present in Fig. [T]three cooling histories for e = 1 
(Wilson action) , e = (slightly overimproved action) , 
and e = — 1 (strongly overimproved action, our choice 
in this paper). The left panel shows this compari- 
son for a gluodynamics configuration taken from the 
confined phase, the right panel for a gluodynamics 
configuration from the transition region. In the full 
statistics analysis of this paper we use e = — I . 

Both action and topological charge were calcu- 
lated with the help of the improved lattice ver- 
sion of the field strength tensor [55J as S/Si nst — 
£x, Q {E a -E- + B-- B a ) /(I6^ 2 ) and Q = £ X>Q E a ■ 
B a /(8ir 2 ). 

As a result, for gluodynamics in the confined phase 
we find calorons to be stable during extended cooling. 
This is well seen in the left part of Fig. [5J One notices 
that the modulus of the topological charge approaches 
some plateau very early. The action (expressed in 
units of instanton action) later converges towards this 
plateau as well. The topological charge plateau is sta- 
ble, what can be understood as stability of calorons, 
which (in the confined phase) are appearing in the 
form of three finite action dyons. 

For gluodynamics in the transition region (early de- 
confined phase), calorons which are rarely present in 
lattice configurations are unstable and cascading down 
to a trivial vacuum as it is seen in the right part of 
Fig. O The first topological charge plateau sometimes 
happens to be different from the first (anti)selfdual 
plateau (this is the case shown in the right panel of 
Fig. [5]). In this temperature range the topological sus- 
ceptibility calculated by the use of topological charges 
measured on the two kinds of plateaus takes different 
values, but the difference is never larger than 6%. 

Deep in the deconfined phase, at T/T c = 1.8, we 
have found that both topological charge and action 
go to zero values without any intermediate plateaus. 

In what follows we will use only topological charge 
values determined on the first topological charge 
plateau. In Fig. [3] again for gluodynamics the Monte 
Carlo time histories of the topological charge are 
shown for the first 200 measurements and the his- 
togram of these values is presented with full statistics. 
These results are presented for the confined phase 
(T = 0.883 T c ) in the upper part and for the decon- 
fined phase (T = 1 .117 T c ) in the lower part of the 
figure. The data indicate that the topological charges 
are well decorrelated. 

We found very similar features of cooling plateaus 
in full QCD. 



IV. TEMPERATURE DEPENDENCE OF 
THE TOPOLOGICAL SUSCEPTIBILITY 

Our results for the topological susceptibility 
f 



x = v (Q 2 ) 



(4) 



for gluodynamics are presented in the left part of Fig. 
IHfor lattices 16 3 ■ 4 (red circles), 24 3 • 4 (blue up tri- 
angles) and 24 3 • 6 (green down triangles). Respective 
results for full QCD are shown in the right part of Fig. 
[U They were obtained by applying cooling to the lat- 
tice ensembles of the DIK collaboration for 16 3 • 8 (red 
up triangles) and for 24 3 • 8 (blue down triangles) . 

From the left part of Fig. |4] we see that the results 
obtained on lattices I6 3 ■ 4 and 24 3 • 4 are close to each 
other. This means that finite volume effects are prac- 
tically absent. The comparison of results from lattices 
16 3 -4 and 24 3 -6 shows that discretization effects in the 
confinement phase are small while sizeable discretiza- 
tion effects are present in the deconfined phase. In 
section [V] we suggest an explanation of this feature. 
On the right panel of Fig. 2]- devoted to full QCD - 
one can see that x is more or less constant below T c 
and starts to decrease slowly above T c . The data also 
indicates finite volume effects to be small, similarly to 
the gluodynamics case. 

The topological susceptibility calculated in this 
work for gluodynamics is compared in Fig. [5] (left) 
with the topological susceptibility calculated by the 
use of the index theorem for a chirally improved Dirac 
operator without cooling [4q |. For full QCD the com- 
parison with results obtained by the use of the in- 
dex theorem for the overlap Dirac operator [47| is 
presented in Fig. [5] (right). In the case of gluody- 
namics we find rather good agreement for all data 
points (uncertainties in the determination of T/T c 
should be taken into account) apart from one point 
at T/T c p» 0.95. The most essential qualitative differ- 
ence between our results and those of Ref. [46| is that 
our data show a sharp drop of x at the transition while 
the data of Ref. [46[ indicate a much smoother behav- 
ior. For full QCD the results agree in the confinement 
phase but disagree in the deconfincment phase. 

In Fig. [5] we compare our gluodynamics results with 
the early results of Ref. [45[ (blue up triangles). In 
this paper the field theoretical method was applied 
for measuring the two-point function of the topolog- 
ical density. This method circumvents the problem 
of defining a topological charge for each configuration 
and uses multiplicative and additive renormalization 
in order to relate the two-point function at zero mo- 
mentum to the continuum topological susceptibility. 
The comparison is made between our results for lat- 
tices 24 3 ■ 6 (green down triangles) and 16 3 ■ 4 (red 
circles) and the results of Ref. [45j which refer to a 
32 3 • 8 lattice (blue up triangles). The agreement with 
our results obtained on finer 24 3 • 6 lattices looks very 
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FIG. 1: Cooling histories of \Q\ corresponding to e = 1 (Wilson action, red dotted line), e — (slightly overimproved 
action, blue dashed line) and e — —1 (overimproved action as used in the rest of this paper, black solid line). Left: different 
cooling histories for one configuration from the confined phase. Right: different cooling histories for a (conditionally 
stable) caloron in a configuration taken from the transition region. 
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FIG. 2: Typical cooling histories of configurations from the confined phase (left). The cooling history of an (finally 
unstable) caloron in a configuration taken from the transition region (right). The red dashed lines represent \Q\, the 
blue solid lines represent S. 



good in both phases. This could be an indication of 
the absence of discretization effects on our L t = 6 
lattices. This is supported by comparison with very 
recent results of Ref. [56| where few sweeps of cool- 
ing were used to evaluate the topological charge on a 
lattice of size 40 3 • 10 . 



V. INTERPRETATION IN TERMS OF A 
DYONIC PICTURE 

Let us interpret our findings in terms of the caloron 
or dyon content of the gluonic field ensembles and of 
the behavior of the average Polyakov loop (related to 
the holonomy of the caloron configurations) during 
the cooling process. 

We discuss the case of gluodynamics first. The ac- 
tion of a caloron does not depend on the holonomy 
(see Appendix and Fig. [5] for an illustration of this 
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FIG. 3: The distributions of topological charges (left) and 
the Monte Carlo time histories along the Monte Carlo 
chain (right) of pure gluodynamics are shown, for the con- 
fined phase in the upper part and for the deconfined phase 
in the lower part. Confinement is represented by a sample 
at T — 0.883 T c , deconfinement is illustrated by a sample 
at T — 1.117 T c , in both cases with an actual statistics 
of 1000 L t = 4,L a — 16 lattice configurations (only partly 
shown in the history). 



fact). So, while the holonomy changes considerably in 
the process of cooling, the action of the (multi)caloron 
(gradually formed in the lattice field configurations by 
cooling in the confined phase) is not changing. Due to 
this stability of the (multi)caloron action the holon- 
omy has no preferred direction to evolve (see the left 
panel of Fig. [7] referring to the confined phase before 
and after cooling). In the deconfined phase on the 
other hand, even if the three center sectors are equally 
represented (see the right panel of Fig. [7] referring to 
the deconfined phase before and after cooling), under 
cooling the holonomy always moves towards the corre- 
sponding corner of the plot, since this is the direction 
where asymmetric dyon-antidyon pairs according to 
the experience in the SU(2) case (see the Introduc- 
tion) could minimize their total action. 

We compare this with full QCD. In this case, cen- 
ter symmetry is slightly violated by the dynami- 
cal fermions already at low temperature ("confined 
phase"). The left panel of Fig. [5] shows where the 
holonomy moves to in the result of cooling. We see 
how the small direct violation of center symmetry in 
the confined phase of full QCD is amplified by cooling 
(again due to the effect of minimization of the action 
of asymmetric dyon-antidyon pairs) . The right panel 
of Fig. [8] shows the holonomies before and after cool- 
ing in the deconfined phase. Here direct and sponta- 
neous violation of central symmetry are summed and 
the asymmetry in the evolution of holonomy is more 
profound. 

The stability of calorons in the confined phase and 
their partial instability in the deconfined phase, both 



in the case of gluodynamics, can be understood if one 
takes into account that during the process of cool- 
ing the holonomy remains nontrivial in the confined 
phase whereas it rapidly becomes trivial in the de- 
confined phase. In the confined phase the caloron, 
for nontrivial holonomy consisting out of three local- 
ized dyons, should first form its nondissociated state 
by recombination of three dyons before it could "drop 
through the meshes" . This, however, is improbable in 
the process of overimproved cooling due to the repul- 
sion between dyons |22j | that the overimproved action 
induces. 

In the deconfined phase of gluodynamics with 
holonomy being almost trivial in the cooling process, 
two light, delocalized dyons forming a caloron to- 
gether with a heavy, localized dyon are already over- 
lapping with the latter. This could be the reason 
of calorons being finally unstable in the deconfined 
phase. 

In QCD, for comparison, the transition from stabil- 
ity to instability develops not so sharply and this is 
connected to the less rapid change with temperature 
of the evolution of the holonomy during cooling as it 
can be seen in Fig. HI 

As we have already seen, the comparison of results 
from lattices 16 3 * 4 and 24 3 * 6 shows the presence of 
sizeable discretization effects in the deconfined phase. 
They can be interpreted as follows. We have discussed 
the partial instability of calorons in the deconfined 
phase. It is quite natural to expect that calorons 
are more stable on a finer lattice where more effort is 
needed to let them "fall through the smaller meshes" 
of the lattice. Hence, on finer lattices the topological 
susceptibility obtained by cooling in the deconfined 
phase should be larger. 



VI. CONCLUSIONS 

In this work we have studied numerically, across 
the phase transition in gluodynamics and across the 
crossover in full QCD, the topological objects, that re- 
main conserved or decay after a suitably chosen cool- 
ing procedure. We have applied the cooling method 
to find the topological content of the first topological 
charge plateau and of the first (anti)selfdual plateau 
in the cooling history of lattice configurations. 

In gluodynamics, we have been able to recognize 
the phase transition point T c as the point separating 
(within a rapid transition) temperature regions with 
and without surviving (anti)selfdual (topological) ob- 
jects. We compare our observations for gluodynamics 
with full QCD with Nf = 2 dynamical flavors in the 
vicinity of the crossover temperature T c . 

We are convinced that the (anti)selfdual plateaus 
can be fully characterized by (multi)caloron configu- 
rations as discussed in [2J|. In the rare cases, in the 
transition region, where the first topological charge 
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plateau could differ from (anti)selfdual plateaus it is 
natural to assume an admixture of caloron-anticaloron 
or dyon-antidyon pairs. 

We calculated the topological susceptibility based 
on the topological charges of individual lattice con- 
figurations identified in this way. We investigated the 
transition or crossover, respectively, from the confined 
to the deconfined phase with the help of this topolog- 
ical susceptibility %. 

In pure gluodynamics with increasing temperature, 
X turned out to drop sharply down from a value close 
to the zero-temperature value just at the (first order) 
phase transition, the latter determined from the be- 



havior of the Polyakov loop or other quantities. Our 
finding tells us that the topological susceptibility de- 
termined via (overimproved) cooling can be used as 
an alternative indicator for the transition itself. 

At the same time the drop off of \ could be viewed 
in terms of stable (multi)calorons on the confine- 
ment side and in terms of unstable dyon-antidyon 
pair or caloron configurations on the deconfincment 
side. The different behavior on both sides could 
be explained as triggered by the holonomy directly 
related to the 3-space averaged Polyakov loop for 
each gauge field configuration. We demonstrated that 
(over)improved cooling drives the holonomy to non- 
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FIG. 7: For gluodynamics, scatter plots of holonomy (spatially averaged Polyakov loop) for thermalized configurations 
and after cooling. Left panels: for the confined phase (T/T c — 0.97) and at the corresponding stable, non-trivial 
(anti)selfduality plateaus; right panels: for the deconfined phase (T/T c — 1.09) and correspondingly after 1500 cooling 
sweeps at (almost) vanishing action values. 
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FIG. 8: For full QCD, scatter plots of holonomy (spatially averaged Polyakov loop) for thermalized configurations and 
after cooling at stable, non-trivial (anti)selfduality plateaus. Left panels: for the confined phase (T/T c — 0.96); right 
panels: for the deconfined phase (T/T c = 1.05). 



trivial (trivial) values in the confinement (deconfine- 
ment) phase with the consequence of getting stable 
(unstable) (multi)caloron configurations due to their 
symmetric (non-symmetric) dyon content. 

In full QCD the drop of the topological suscep- 
tibility turned out to be softer and to set in at a 
temperature that is slightly above the pseudocritical 
temperature T c determined by the maximum of the 
Polyakov loop susceptibility [49| . Also in this case the 
former discussion of stable (below the crossover) or 
unstable (above the crossover) (multi)caloron config- 
urations applies, since the holonomy behaves similarly 
in spite of the Z(3) breaking effect by the dynamical 
fcrmion degrees of freedom taken into account. 

We compared our results on topological susceptibil- 
ities with those of other authors employing different 
methods and have found a reasonable agreement. 

All together, our results give us some confidence 
that large-scale topological objects play a major role 
in the change of thermal gauge field ensembles at the 
deconfincmcnt phase transition in gluodynamics as 
well as at the crossover phenomenon in full QCD. 
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Appendix: 5(7(3) calorons 

The SU(N) instantons at finite temperature (or 
calorons) with non-trivial holonomy [ 1 ll — 1 1 3j | can be 

considered as composites of N constituent monopoles, 
seen only when the Polyakov loop at spatial infin- 
ity (holonomy) is non-trivial. In the periodic gauge, 
A fJ ,(t + f3,x)=A^(t,x) it is defined as 



Voo = lim P exp( / A 

l J h°° Jo 



{x,t)dt). 



(5) 



After a suitable constant gauge transformation, the 
Polyakov loop can be characterised by real numbers 
A i m=i,...,n Em=iMm = 0) that describe the eigenval- 



ues of the holonomy 

P%o = exp[27ridiag(/ii,...,/i„)], 
Hi < ... < n n < /i n +i = l + A«i- 



(G) 



In units, where the inverse temperature f3 = 1, a sim- 
ple formula for the SU(N) action density can be writ- 
ten 0, [H : 

TrF*,(x) = %%hgil>(x), (7) 

ip(x) = 7;tr{A n • • -Ai) - cos(27rt), 



a _ J_ fr m \vm-y m +r\ fc m s, 
m ~ r m \ r m+1 J [s m c r 



with r m = \x— y m \ and y m being the center of mass radii 
of m constituent monopoles, which can be assigned a 
mass 87r 2 i/ m , where v m = Hm+i — /i m . Furthermore, 
Cm = cosh(27ri/ m r m ), s m = sinh(27w m r m ), r n+1 = i\ 
and y„+ i=yi. 

For SU{2>) calorons we correspondingly para- 
metrize the asymptotic holonomy as V^ = 
diag(e 27rw ,e 27rV2 ,e 27rV3 ), with fix < \i 2 < H3 < fi± = 
1 + fix and fix + A«2 + A*3 = 0. Let j/i, yb and i/3 be 
three 31? position vectors of dyons remote from each 
other. Then a caloron consists of three lumps carrying 
the instanton action split into fractions mi = /i 2 — fJ>\, 
m 2 = H3 ~ A*2 and 7713 = ^4-/13, concentrated near 
the £. 



Provided the constituents are well separated, the 
Polyakov loop values at their positions y mi m = 
1,2,3 are [57| 



V(yi) = diag( 



V{y 2 ) 



diag( 
diag(^ 



27Tl/Ji 



^Trifj.2 



e 27rV3 ), 
e' mtil ), (8) 



The complex numbers representing the trace of 
Polyakov loop PL = iTr"P occupy some region on 
the complex plane (see e.g. Fig. [7]). The holon- 
omy V^o is traced to point PL^ that is close to 
(one third of) the trace of Polyakov loop averaged 
over all lattice points and is near zero in the con- 
fining phase and near one in the deconfining phase. 
There is a one-to-one correspondence between (one 
third of) the trace of Polyakov loop and three numbers 
mi, m2, m 3j TUi + m 2 + 7713 = 1 that define the eigen- 
values of SU(3) matrix and hence its trace. Three 
numbers can be represented by the inner point of reg- 
ular triangle for which the sum of the lengths of three 
perpendiculars to triangle sides is constant (equal to 
one, see Fig. |9|). The region on the complex plane 
occupied by (one third of) the trace of Polyakov loop 
can be considered as some nonlinear deformation of 
this regular triangle. The point O on Fig. [9] corre- 
sponds to PLqo while points Ai,A 2 , A 3 correspond to 
the values of (one third of) the trace of the Polyakov 
loop at the constituent positions ©• 
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